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Abstract—Bearing voltages and corresponding currents in 
electric machines driven by pulse width modulation (PWM) 
converters with fast switching and high dv/dt can cause prema- 
ture bearing failures. This paper evaluates the bearing voltage 
reduction by using insulated shafts and bearings. An equivalent 
circuit representation of electric machines taking into account 
high-frequency effects is developed to show the production 
mechanism of bearing voltage, based on which simulations are 
performed with detailed finite element models for transient 
and simplified equivalent circuit for steady-state analysis. The 
steady-state equivalent circuit is then calibrated, by combined 
numerical calculations and experimental measurements, and used 
to predict the bearing voltage for various scenarios, proving the 
effectiveness of using insulated shaft and bearings in reducing 
the steady-state bearing voltage. 

Index Terms—Bearing current, bearing voltage, capacitance, 
electric machine, insulated bearing, insulated shaft, PWM con- 
verter. 


I. INTRODUCTION 


Bearing failures in electric machines and electric machine- 
driven equipment are attributed to mechanical/thermal causes 
or bearing currents. Shaft voltages and the resulting bear- 
ing currents in line-operated electric machines are mainly 
produced by the homopolar flux, magnetic asymmetry and 
electrostatic effect. During the past decades, with the prolif- 
eration of variable speed drives using pulse width modulation 
(PWM) voltage source inverters, a growing number of bearing 
current induced premature bearing failures have been reported, 
calling for more research effort into bearing currents and the 
associated mitigation techniques [1]. 

All rotating machines potentially have bearing current is- 
sues, either of large or small horsepowers [2]. It has been 
shown that small inverter-fed ac machines are likely to suffer 
from electrical discharge machining (EDM) bearing currents, 
and larger machines are likely to be exposed to high-frequency 
circulating bearing currents [3]. 

The voltage across the bearing lubricant film, i.e., the 
bearing voltage Vp, can be built up by the electrostatic charge 
accumulation, the common mode voltage (CMV) through 
capacitve coupling, and the common mode current or magnetic 
asymmetry through inductive coupling. The bearing current 
will flow and damage the surface of bearing raceways once 


the instantaneous bearing voltage exceeds a given threshold 
value [2]. 

Main bearing current types include: the localized current 
induced by homopolar fringing magnetic fluxes, the leakage 
current to the ground, the non-circulating type caused by 
the shaft-to-frame voltage and the circulating type introduced 
by the shaft end-to-end voltage [4]. The bearing voltage 
is typically used as an indicator of bearing failures since 
both the non-circulating and circulating bearing currents are 
closely related to the potential drop across the bearing [5], [6]. 
The accurate prediction of bearing voltage and cost-effective 
mitigation solutions are important topics in this research area, 
such as the recent work [7], [8]. 


II. EQUIVALENT CIRCUIT REPRESENTATION OF ELECTRIC 
MACHINES FOR BEARING VOLTAGE STUDY 


This paper takes a 3-phase 12-slot 10-pole permanent- 
magnet synchronous machine (PMSM) as an example study 
for the bearing voltage and evaluates its reduction potential 
by employing insulated shaft and bearings based on cali- 
brated models, attempting to develop an effective and efficient 
modeling approach for bearing voltage and current study. By 
”insulated shaft”, the author refers to the shaft whose diameter 
is reduced to fit in an insulation layer, e.g., plastic, between 
the shaft and inner ring of the bearing. Insulated bearing” 
refers to the bearing whose outer surface of the outer ring 
and/or inner surface of the inner ring is E-coated or covered 
by plastic. 

The equivalent circuit representation of the PMSM under 
study is developed for bearing voltage simulation, considering 
both the low-frequency and high-frequency behaviors, as illus- 
trated in Fig. 1. The bearing voltage waveform shown in Fig. 
1 depends heavily on the inverter output, cable impedance, 
ground impedance, high-frequency parameters of the machine 
and bearing parameters. 

In order to show the effect of distributed parameters in 
transient simulations, the equivalent circuit includes details on 
winding turns. The construction of the 12s10p PMSM illus- 
trated in Fig. 2 has an outer rotor with surface-mounted bonded 
ferrite magnets and an inner stator with 12 concentrated tooth 
wound coils, as shown in Fig. 3. 
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Fig. 1. Equivalent circuit for the PMSM under test taking into account both the low-frequency behavior and high-frequency effects. Windings are represented 
turn by turn and each winding turn block includes a resistance, an inductance and a back electromotive force (EMF). Only capacitances related to Phase A 


coils are shown for clear illustration. 
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Fig. 2. Exploded view of a PMSM with an outer rotor. 
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(a) 
Fig. 3. PMSM under test: (a) outer rotor, (b) inner stator. 


There are 39 turns in each coil and two 0.352mm magnet 
wires in hand per turn. Both full ceramic bearings and steel 
ball bearings are used in the test for model calibration. Photos 
of the full ceramic bearing 608 and steel ball bearing 608 are 
included in Fig. 4. 

In electric machines, the low-frequency behavior involves 
resistances and inductances at the fundamental frequency, 
which are commonly used for motor operation simulation 
and control purpose. The bearing voltage shown in Fig. 1, 
by contrast, is mainly induced by distributed capacitances 
and closely related to high-frequency inductances, resistances 
and parasitic capacitances of windings. These hig-frequency 
parameters need to be calculated and measured in ways 


(a) 


(b) 
Fig. 4. Full ceramic bearing (a) and steel ball bearing (b). 


different from those for resistances and inductances at the 
fundamental frequency. 


III. BEARING VOLTAGE SIMULATION 
A. Problem Formulation 


Existing methods for bearing voltage prediction are mainly 
based on equivalent circuits/networks, where the circuit pa- 
rameters need to be obtained first. The dc resistance is usually 
calculated from the wire gauge and average turn length. The 
influence of temperature can also be taken into account by us- 
ing the conductivities at different operating temperature. FEA 
has been widely used in estimating machine parameters at the 
fundamental frequency, such as d-axis and q-axis inductances 
for PM machines [9]. 

The measurements of high-frequency winding resistance, 
inductance and parasitic capacitance are typically performed 
by impedance analyzers or the equivalent in a wide frequency 
range. Studies have also been reported on the extraction of 
electric machine parameters at high frequencies using FEA 
[10]-[13]. Considering the fact that high-frequency resistances 
and inductances in electric machines are frequency-dependent, 
methods are required that can take into account the frequency 
dependence of parameters. 


B. FE Modeling 


This paper proposes to simulate the transient response of 
bearing voltage in electromagnetic FE transient solver with 


Fig. 5. Detailed solid model for the FE modeling of the motor under study 
for parameter extraction and transient analysis. Stator windings are modeled 
turn by turn to capture the distribution effect of parameters. 
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Fig. 6. Phase resistance and inductance vs. frequency from FEA. Eddy 
currents in conductors are modeled whereas the coreloss effect on the field is 
not considered. 


the help of external circuits. The distributed resistance and 
inductance are modeled coil by coil. Capacitances need to be 
extracted first and then incorporated into external circuits. To 
capture the bearing voltage transients, machine windings need 
to be modeled in detail to consider the distribution effect of 
conductors. Each in-slot conductor or turn can be treated as 
an independent node to consider the voltage distribution along 
the wire [14]. 

High-frequency resistances and inductances can be calcu- 
lated explicitly by an eddy current solver or implicitly by a 
magnetic transient solver. The FE model with winding details 
of the PMSM under study is shown in Fig. 5. The calculated 
phase resistance and inductance are frequency-dependent, as 
illustrated in Fig. 6. In the range where the frequency is 
lower than 100kHz, the resistance and inductance are nearly 
constant. With futher increase in the frequency, the resistance 
increases drastically and the inductance reduces gradually due 
to the eddy effect in conductors. 

Capacitances in electric machines can be calculated using 
an electrostatic solver. The 3D mesh for the electrostatic model 
for the calculation of capacitances shown in Fig. 7 has a 
total 29.14 million tetrahedrals. The calculated capacitances 
are independent of the frequency and the voltage distribution 
is visualized in Fig. 8. 


Fig. 7. Meshing of the 3D FE model for an electrostatic solver in ANSYS 
Maxwell [15]. There are in total 29.14 million tetrahedral elements and it 
took 35 hours to solve the problem on a high-performance PC with an 8-core 
Intel Xeon processor and 128GB RAM. 


Fig. 8. Voltage distribution within the motor when 1V is applied to the stator 
winding and the stator core is grounded. 


In order to show the capacitances in machines with winding 
details, a 2D FEA model is built with the conductor layout 
shown in Fig. 9. Coils are assumed to be wound regularly with 
a needle winder starting from the bottom of the slot. There are 
3 conductor layers and the two wires in hand are placed in 
pairs. The capacitance matrix shown in Fig. 10 has 4 traces 
highlighted by light colors. The two light blue traces denote 
capacitances between adjacent conductors from the same layer. 
The two yellow traces denote capacitances between conductors 
from adjacent layers. The influence of conductor positions on 
capacitances can also be observed from results shown in Fig. 
11. 

For steady-state bearing voltage analysis, all the winding 
turns are equipotential and the equivalent circuit shown in Fig. 
1 can be simplified to a pure capacitve network illustrated in 
Fig. 12. The constructions of stator and rotor in Fig. 3 show 
that the laminated stator core and the bearing holder are part 
of the same assembly and can be treated as a single node. 
Similarly, the rotor back iron, PM ring and the shaft can be 
regarded as another node. The winding represents the third 
node since it is insulated from the stator core and rotor. 


The steady-state bearing voltage mirrors the CMV and is 


Fig. 9. Layout of turns in one coil of a phase winding. There are two 0.3mm 
wires in hand and coils are assumed to be wound regularly with a needle 
winder for the simplicity of modeling and analysis. 
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Fig. 10. Capacitance matrix for one coil of Phase B winding from 2D FEA. 


2 10 
- E 
= wn 
515 ae £ 
z 0.6 2 
oO oO 
o 1 © 
8 0.4 8 
O £ 
g 0.5 028 
(dp) 1 
ó ð 0.0 S 
0 5 10 15 20 25 30 35 40 
Turn No. 


Fig. 11. Winding turn related capacitances from 2D FEA. 


determined by the voltage divider, which can be written as: 
Cwr n 

Cwr + Csr + (Cor + Coz) 

where Cwr is the lumped capacitance between stator wind- 


ing and rotor, Csr the lumped capacitance between stator core 
and rotor. Cy; and Cyz are the capacitances of bearing 1 and 


V, = 


CMV, (1) 


Node 
(Stator core + bearing holder) 


Node 
(winding) 


Fig. 12. Equivalent circuit for steady-state analysis. 


TABLE I 
COMPARISON OF CALCULATED AND MEASURED CAPACITANCES IN THE 
MOTOR WITH FULL CERAMIC BEARINGS MADE OF ZrO2. THE 
DIELECTRIC CONSTANT OF Z,;O2 HAS A SMALL IMPACT ON THE FE 
PREDICTED CAPACITANCES. 


Individual FEA FEA Meas. Percentage 
capacitors €r =10 [pF] er =23 [pF] [pF] Difference [%] 
Cws 316.3 316.2 286.2 +10.5 
Csr 64.9 66.9 65.1 -0.3/+2.8 
Cwr 4.5 4.5 4.0 +12.5 


bearing 2, respectively. It can be seen that the steady-state 
bearing voltage is independent of the capacitance between 
stator winding and stator core Cws, as well as the capacitances 
between phase windings Cy. 


IV. DETERMINATION OF BEARING CAPACITANCES BASED 
ON A COMBINED NUMERICAL AND EXPERIMENTAL 
TECHNIQUE 


The winding-to-stator-core capacitance Cws, stator-core-to- 
rotor-core capacitance Csr and winding-to-rotor-core capac- 
itance Cwr are calculated directly from 3D FEA shown in 
5, including the bearings. Full ceramic bearings are used in 
the FE model to reduce the computation time. The calculated 
and measured values are tabulated in TABLE I, showing a 
maximum difference of 12.5%. 

The capacitance of a steel ball bearing depends on its 
construction, lubricant property, lubrication status, operating 
temperature, speed and load. The dielectric constant of the 
lubricant mixture (ep), lubrication status (€e) and additional 
capacitance (Cother) introduced by metal ball retainers and 
seals are difficult to be determined without experimental 
measurements. To obtain accurate capacitance for steel ball 
bearings, a combined numerical and experimental approach is 
developed, which models the steel ball bearing capacitance 
based on analytical equations according to the Hertzian con- 
tact theory and uses measured capacitances to determine the 
unknown coefficients. 

Experimental tests were conducted on the PMSM under 
study at 1,200r/min with a fixed DC bus voltage of 24V. 
The space vector PWM (SVPWM) based field-oriented control 
was employed to control the motor speed with a switching 
frequency of 20kHz. The CMV and bearing voltage waveforms 
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Fig. 13. Measured common mode voltage and bearing voltage at 1,200r/min 
for: (a) full ceramic bearing, (b) steel ball bearing. 


were recorded by an oscilloscope with a sampling frequency 
of 50MHz, as shown in Fig. 13. The peak-to-peak values of the 
CMV and bearing voltage were measured from the waveforms 
and used to calculate the combination C’,,.+Cy1+Cp2, as shown 
in TABLE II based on the circuit illustrated in Fig. 12. 

It can be seen that the bearing voltage mirrors the CMV 
very well in the tested speed range from 1,200r/min to 
3,360r/min. Voltage ringing is observed at the rising and falling 
edges of the bearing voltage, which might be introduced by 
the resistances and inductances neglected in the steady-state 
analysis. The dc bias in the CMV caused by the dc bus voltage 
pulsation is eliminated in the bearing voltage. The calculated 
capacitances using analytical equations, FEA and the values 
calibrated by experimental results are also included in TABLE 
II, showing satisfactory agreement. 


V. BEARING VOLTAGE REDUCTION BY INSULATED 
SHAFTS AND BEARINGS 


The potential of bearing voltage reduction by insulating 
the bearing from the shaft and/or bearing holder was studied 
based on the calibrated equivalent circuit model. The dielectric 
constant of the E-coating layer and plastic are assumed at the 
typical value of 3.0. 

A layer of insulation was inserted between the inner ring 
of bearing 608 and shaft or the outer ring of bearing 608 and 
bearing holder. The influence of the insulation thickness on the 
voltages of bearing 608 and 688 shown in Fig. 14 indicates 
that with the increase of the insulation thickness, the voltage of 
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Fig. 14. Bearing voltage reduction versus the thickness of insulation between: 
(a) the inner ring of bearing 608 and shaft, (b) the outer ring of bearing 608 
and bearing holder. Inner rings have smaller surfaces and therefore produces 
larger additional impedance and voltage reduction. 


the insulated bearing was reduced but the voltage of the non- 
insulated bearing increased slightly. An 80% and 55% voltage 
reduction can be achieved by insulating the inner ring and 
outer ring, respectively. Insulating bearing 688 shows similar 
trends. 

Fig. 15 shows that both bearing voltages can be reduced by 
insulating bearing 608 and bearing 688 at the same time. With 
a 1-mm thick insulation layer inserted between the inner rings 
and shaft, the bearing boltages can be reduced by around 70%. 
In addition, insulating the inner rings reduces more bearing 
voltage due to the smaller facing surfaces. 


VI. CONCLUSION 


This paper studies the bearing voltage reduction potential 
of using insulated shaft and bearings based on the equivalent 
circuit representation and FEA. To improve the accuracy of 
bearing voltage prediction, a model calibration process by 
combining numerical calculations and experimental measure- 
ments is developed. 

In order to track the dynamics in bearing voltages when 
the motor is running, FE models with detailed winding infor- 
mation are required. For steady-state analysis, the simplified 
equivalent circuit may be used. With the combined numerical 
and experimental approach, the percentage error for bearing 
capacitance calculations can be reduced to less than 5%. 
Insulating a single bearing can reduce the voltage on the 


TABLE II 
COMPARISON OF CALCULATED AND MEASURED CAPACITANCES. THE REGRESSION ANALYSIS SHOWS THAT THE OPTIMAL COEFFICIENTS ARE €f = 
5.0, €e = 1.0, Cother = 0. FE-PREDICTED CAPACITANCE C'sr IS 61.8 PF. 


Full ceramic bearing 


Steel ball bearing 


Speed [r/min] 1,200 1,200 2,400 3,360 
Meas. CMV [V] 24.0 24.0 24.2 24.4 
Meas. V, [V] 1.80 1.14 1.22 1.34 
Meas. Csr + Coi + Coz [pF] 65.1 105.15 98.0 88.85 
Calc. Cy1 [pF] 2.40 (FE) 21.2 (anal.) 16.3 (anal.) 14.6 (anal.) 
Calc. Cy2 [pF] 2.29 (FE) 21.5 (anal.) 16.4 (anal.) 11.3 (anal.) 
Calc. Csr + Coi + Co2 [pF] 66.5 104.5 94.5 86.9 
Percentage difference [%] 2.1 - 0.6 3.6 22: 
REFERENCES 
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Fig. 15. bearing voltage reduction versus thickness of insulation layers on: (a) 
inner surfaces of both bearing 608 and 688, (b) outer surfaces of both bearing 
608 and 688. Voltages of bearing 608 and 688 can be reduced effectively by 
insulating the inner or outer surfaces of both bearings simultaneously. 


insulated bearing but will adversely increase the voltage on the 
non-insulated bearing. To reduce the voltage of both bearings, 
it is recommended to insulate the inner rings of both bearings. 
A maximum voltage reduction of 70% is predicted in the 
PMSM under study. Studies for other machine ratings are 
underway. 
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